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MATCHING PARTON SHOWER AND MATRIX ELEMENTS IN QED 
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We report on a high-precision calculation of the Bhabha process in QED, of interest for precise 
luminosity determination of low-energy electron-positron colliders. The calculation is based on 
the matching of exact next-to-leading order corrections with a Parton Shower algorithm. The 
structure of the algorithm (implemented in an improved version of the event generator BABAYAGA) 
is illustrated, with a discussion on the resulting theoretical uncertainty, of the order of 0.1%. 
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1. Introduction 

The measurement of the ratio R 



^T(e^ 



hadrons)/(T(e+e 



) at 



flavour factories, such as DA$NE, VEPP- 
2M, BES, KEK-B, PEP-II and CLEO, is of 
primary importance for the precise determi- 
nation of the anomalous magnetic moment 
of the muon, (17 — 2)^, and of the running 
of the electromagnetic coupling aQEoiQ^)- 
The cross section values entering R are af- 
fected by the uncertainty on the knowledge 
of the machine luminosity, which is, in turn, 
related to the uncertainty on the theoretical 
knowledge of the cross section of a reference 
QED process, typically large angle Bhabha 
scattering, jX^ and 77 production. To 
keep under control such an uncertainty, high- 
precision calculations of these QED processes 
and relative Monte Carlo generators, are re- 
quired. Large-angle Bhabha scattering is of 



particular interest because of its large cross 
section and its clean experimental signature. 
To simulate the experimentally relevant dis- 
tributions and calculate the cross section of 
the Bhabha process, KLOE and CLEO col- 
laborations make use of the QED Parton 
Shower (PS) generator BABAYAGA, developed 
in Refs. with a precision of 0.5%. At 
present a reduction of such a theoretical sys- 
tematics is needed for several reasons. First, 
the experimental luminosity error quoted by 
KLOE is presently 0.3% ^. Secondly, the 
measurement of the hadronic cross section in 
the 7r+7r~ channel at VEPP-2M has achieved 
a total systematic error of 0.6—1% in the re- 
gion of the p resonance ^, which requires, in 
turn, an assessment of the collider luminos- 
ity at the level of 0.1%. Last but not least. 
Charm and _B-factories will reach in the near 
future the error of 1% in the luminosity mea- 
surement. 
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At the 0.1% level the non logarith- 
mic contributions present in exact next-to- 
leading (NLO) pcrturbativc calculations as 
well as the resummed leading logarithmic 
contributions taken into account in the PS 
approach are expected to be relevant. A 
matching algorithm which allows to incor- 
porate the next-to-leading terms within the 
PS structure of an event generator such as 
BABAYAGA, without double counting at first 
order in a of the leading corrections already 
accounted for by the PS, has been developed 
in Ref. ^ and will be reviewed in the follow- 
ing. An estimate of the remaining theoretical 
uncertainty is also discussed. 

2. Matching NLO corrections 
with Parton Shower 

A general expression for the cross section 
with the emission of an arbitrary number of 
photons, in leading-log (LL) approximation, 
can be cast in the following form: 

oo 

da21 = n(Q2,e) V- \Mn,LL\^ (1) 

' n! 

where Il{Q^,e) is the Sudakov form-factor 
accounting for the soft-photon (up to an en- 
ergy equal to s in units of the incoming 
fermion energy E) and virtual emissions, e 
is an infrared separator dividing soft and 
hard radiation and is related to the en- 
ergy scale of the process. |A^n,LLp is the 
squared amplitude in LL approximation de- 
scribing the process with the emission of n 
hard photons, with energy larger than e in 
units of E. d^n is the exact phase-space ele- 
ment of the process (divided by the incoming 
flux factor) , with the emission of n additional 
photons with respect to the Born-like final- 
state configuration. 

According to the factorization theorems 
of soft and/or coUinear singularities, the 
squared amplitudes in LL approximation can 
be written in a factorized form. In the fol- 
lowing, for the sake of clarity and without 



loss of generality, we write photon emission 
formulas as if only one external fermion radi- 
ates, being the generalization to the real case 
straightforward when including the suited 
combinatorial factors. The one-photon emis- 
sion squared amplitude in LL approximation 
can be written in factorized form as 

\Mi,LL ^ = —- I{k) Mo ^ T;2-Ti V 

(2) 

where 1 — 2 is the fraction of the fermion en- 
ergy E carried by the photon, k is the photon 
four-momentum, I{k) is a function describ- 
ing the angular spectrum of the photon. 

The cross section calculated in Eq. (1) 
has the advantage that the photonic correc- 
tions, in LL approximation, arc resummed up 
to all orders of perturbation theory. On the 
other side, the weak point of the formula (1) 
is that its expansion at 0{a) does not coin- 
cide with an exact 0{a) (NLO) result, being 
its LL approximation. In fact 

dalj^ = [l + C„,ll] |Mo|'d$o + IMi.LLpd^i 

(3) 

whereas an exact NLO cross section can be 
always cast in the form 

da« = [1 + C„] |Mop(i*o + (4) 

The coefficient contains the complete vir- 
tual 0{a) and the 0{a) soft-bremsstrahlung 
squared matrix elements, in units of the Born 
squared amplitude, and |A^i|^ is the exact 
squared matrix element with the emission of 
one hard photon. We remark that Ca.LL 
has the same logarithmic structure as Ca 
and that llP has the same singular be- 
haviour of |A^ip. 

In order to match the LL and NLO calcu- 
lations, we introduce the correction factors, 
which are by construction infrared safe and 
free of collinear logarithms, 

Fsv = 1 + (Ca - Ca,LL) , 

= l + (5) 
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and we notice that the exact 0{a) cross sec- 
tion can be expressed, up to terms of 0{q.^), 
in terms of its LL approximation as 

rf<7a = Fsvil + Co,,LL)\Mo\''d^o 

+ FH\Mi,LL\^d<i>i (6) 

Driven by Eq. (6), Eq. (1) can be improved 
by writing the resummed cross section as 

matched = ^SV 11(0^, e) 
n=0 ■ \i=0 / 

The expansion at 0{a) of Eq. (7) coincides 

now with the exact NLO cross section Eq. (4) 
and all higher order LL contributions are the 
same as in Eq. (1). 

Eq. (7) is our master formula for the 
matching between the exact 0{a) calciilation 
and the QED resummed PS cross section, 
according to which we also generate events. 
The correction factors of Eq. (5) can in prin- 
ciple make the differential cross section of 
Eq. (7) negative in some point, namely where 
the PS approximation is less accurate. Nev- 
ertheless, we verified that this never happens 
when considering typical event selection cri- 
teria for luminosity at flavour factories. 

2.1. Vacuum polarization 

Besides the photonic radiative corrections 
considered above, also the vacuum polariza- 
tion eff'ects must be included in the master 
formula (7), in order to reach the required 
theoretical accuracy. They are accounted for 
by replacing the fine structure constant a = 
a{0) with a{q^) = a/{l — Aa{q'^)), according 
to the algorithm described in Ref. ^. Aa{q'^) 
is the fermionic contribution to the photon 
self-energy: the leptonic and top-quark one- 
loop contributions can be calculated ana- 
lytically in perturbation theory, while the 
remaining five quarks (hadronic) contribu- 
tion, Aa^^^^j,, has to be extracted from data. 
To evaluate Aa\^^^ we use the HADR5N rou- 
tine 6.7. 



3 

Going beyond the Born-like approxima- 
tion, the cross section corrected at 0{a) in- 
cluding also vacuum polarization can be writ- 
ten as ayp = ao^vp + (^sv + '^h^ where agy 
and (T^ are the soft plus virtual and the hard 
photon 0{a) corrections of photonic origin. 
We can go a step further and include vac- 
uum polarization in those terms, in order to 
include also part of the 0{a^) factorizable 
corrections. 

Furthermore, we add to the Born am- 
plitude also the Z exchange diagrams: their 
effect is really tiny and negligible at low en- 
ergies, but can become more important (up 
to 0.1%) around 10 GeV. 

3. Theoretical uncertainty 

Since different implementations of radiative 
corrections beyond exact 0{a) contributions 
differ by higher order effects, a hint of the 
missing radiative corrections which domi- 
nate the theoretical accuracy can be given 
by comparing the predictions of BABAYAGA 
with independent event generators, such as 
BHWIDE^, LABSPV^ and MCGPJ^°. As shown in 
Ref. ^, the results of these comparisons are 
very satisfactory, with differences between 
BABAYAGA and BHWIDE below 0.1% on cross 
sections and ranging up to 1% only in the 
tails of some distributions where the statis- 
tics is low. 

A firmer estimate of the theoretical ac- 
curacy could be given by comparison with a 
complete two-loop calculation, which is how- 
ever not available yet. However, recently 
there has been important progress towards 
the full 0(0^) calculation. At present, two 
different partial contributions have been cal- 
culated: the complete virtual two- loop pho- 
tonic corrections (in the limit » m^, 
with = s, t, u) plus real radiation in 
soft approximation and the virtual Nf = 
1 fermionic contribution inclusive of finite 
mass terms In Ref. ^ the terms in 

BABAYAGA corresponding to these two ap- 
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proximations of the complete O(a^) calcu- 
lation have been extracted and compared, 
showing excellent agreement. The relative 
differences don't exceed the 0.03% level. A 
careful inspection of the analytical expres- 
sions of the differences shows that all log- 
arithmic terms of infrared origin present in 
BABAYAGA have the same coefficients as in the 
two- loop perturbative calculations, with the 
exception of small terms suppressed by pow- 
ers of rrig/s. 

Other O(a^) contributions not consid- 
ered in BABAYAGA (and therefore sources of 
theoretical uncertainty) are the light pair 
corrections and the soft plus virtual 0{a) 
corrections to the real hard radiation. The 
impact of the former contribution has been 
estimated for a sample of typical event selec- 
tion and energies, with f-channel virtual LL 
approximated formulae^'^ and real pair emis- 
sion in soft approximation^^. The impact of 
such corrections has been found below the 
0.05% level. Exact perturbative results for 
the soft plus virtual corrections to real hard 
radiation are not available yet for the com- 
plete s + t Bhabha process, which is of inter- 
est for the case of low energy e+e^ colliders. 
They have been calculated separately for the 
t- and s-channels and their impact studied 
at LEP conditions. Taking into account of 
such experience and that BABAYAGA contains 
all the infrared enhanced terms, the size of 
these 0{a'^) corrections can be estimated to 
be smaller than 0.05%. Adding in quadra- 
ture the perturbative sources of theoretical 
uncertainty would give a theoretical error of 
the order of 0.1% for all considered event 
selections and energies, from 1 to 10 GeV. 
However this value underestimates the to- 
tal error associated with BABAYAGA at Charm 
and S-factories because the routine HADR5N, 
which gives the vacuum polarisation correc- 
tions, produces large errors around the J/ip 
and T resonances, increasing the total theo- 
retical error of BABAYAGA at the 0.2% level in 



those energy ranges. 

F.P. would like to thank the conveners of 
the electroweak parallel session for their kind 
invitation. 
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